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Abstract. We present first results on the quantitative spectroscopic analysis of the photospheric-phase of Type 
II supernovae (SN). The analyses are based on the model atmosphere code, CMFGEN, of Hillier & Miller (1998) 
which solves the radiative transfer and statistical equilibrium equations in expanding outflows under the constraint 
of radiative equilibrium. A key asset of CMFGEN is its thorough treatment of line-blanketing due to metal species. 
From its applicability to hot star environments, the main modifications to the source code were to allow a linear 
velocity law, a power-law density distribution, an adaptive grid to handle the steep H recombination/ionization 
front occurring in some SN models, and a routine to compute the gray temperature structure in the presence 
of large velocities. In this first paper we demonstrate the ability of CMFGEN to reproduce, with a high level of 
accuracy, the UV and optical observations of a sample of well observed Type II SN, i.e. SN1987A and SN1999em, 
at representative stages of their photospheric evolution. Two principal stages of SN are modeled — that where 
iy-^ ' hydrogen is fully ionized, and that in which H is only partially ionized. For models with an effective temperature 

f*"*) , below ~8,000K, hydrogen recombines and gives rise to a steep ionization front. 

The effect of varying the location of the outer grid radius on the spectral energy distribution (SED) is investigated. 
We find that going to 5-6 times the optically-thick base radius is optimal, since above that, the model becomes 
prohibitively large, while below this, significant differences appear because of reduced line-blanketing (which 
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persists even far above the photosphere) and the truncation of line-formation regions. To constrain the metallicity 
and the reddening of SN, the UV spectral region of early-time spectra is essential. We find that the density of the 
photosphere and effect of line blanketing decline as the spatial scale of the SN increases. The density distribution 
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is found to have a strong impact on the overall flux distribution as well as line profiles. For a given base density, the 
faster the density drops, the higher the effective temperature of the model. We also find in cool models that the set 
of Call lines, near 8500A is strongly sensitive to the density gradient. They show a weaker and narrower profile 
for steeper density distributions. Hydrogen Balmer lines are very well reproduced in fully or partially ionized 
models, but underestimated when hydrogen recombines. A reduced turbulent velocity or a flatter density layout 
are found to partially, but not fully, cure this persistent problem in studies of Type II SN. He I lines observed 
in early-time spectra are very well reproduced, even for very modest helium enrichments, likely resulting from 
treatment of important non-LTE effects. At similar early epochs CMFGEN predicts, unambiguously, the presence 
of Nil lines in the blue-wing of both H/3 and He I 5875A. These lines have been observed but so far have generally 
been associated with peculiar emission, from locations far above the photosphere, in the strong adjacent lines. 
Finally, we present a pedagogical investigation on P-Cygni profile formation in Type II SN. Ha is found to form 
very close to the photosphere and thus presents a significant flux-deficit in the red, made greater by the rapidly 
declining density distribution. This provides a clear explanation for the noticeable blue-shift of P-Cygni profiles 
observed in early-time spectra of Type II SN. 

Future studies based on CMFGEN modeling will focus on using Type II SN for the calibration of distances in the 
Universe, as well as on detailed spectroscopic analyses for the determination of progenitor properties. 

Key words, radiative transfer - Methods: numerical - stars: atmospheres - stars: supernovae - line: formation 



1. Introduction 



Send offprint requests to: Luc Dessart, 
e-mail: lucOas . ar izona . edu 



Following the development of radiative transfer algorithms 
suitable for expanding outflows (Mihalas et al. 1975), and 
advances in computer technology, increasingly detailed 
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and successful quantitative spectroscopic analyses of mas- 
sive stars and their winds have been carried out (c.f., 
Schmutz 1997; Hillier & Miller 1998, 1999; Hamann & 
Koesterke 1998). Such quantitative spectroscopic studies 
can shed light on numerous properties of optically-thick 
outflows, with implications for understanding the parent 
object. One can determine the chemical composition, ra- 
dius, luminosity, density, ionization structure and kine- 
matic properties of the star and its wind, as well as the 
temperature of the equivalent blackbody at the thermali- 
sation depth. 

Although stellar winds are radiatively-driven, the 
mechanisms controlling the escape of photons from the 
optically-thick flow are similar in nature to that occur- 
ring in the freely-expanding fast ejecta associated with 
supernovae (SN) explosions, occurring either via core- 
collapse of a massive star (H-free Type Ib/c or H-rich 
Type II; Arnett et al. 1989) or thermonuclear runaway of 
a Chandrasekhar mass white-dwarf (Type la; Hillebrandt 
& Niemeyer 2000). In the context of SN, spectroscopic 
modeling allows inferences on the properties of the pro- 
genitor of the explosion to be made, providing clues to 
the evolution and death of stars, as well as the formation 
of compact objects. 

Additionally, massive stars, and SN explosions to a 
much larger extent, are so intrinsically luminous that at- 
tempts to use them as distance calibrators have been 
made. For example, massive stars can be used to constrain 
the Hubble constant via the Wind-Momentum-Luminosity 
relation (Kudritzki et al. 1999). Light curves of Type la 
SN have been used to claim that we live in an acceler- 
ating universe (Riess et al. 1998). Finally, although Type 
II SN are not as bright as Type la events, they can be 
used with variations on the Baade's method (Baade 1926), 
i.e. the Expanding Photosphere Method (EPM; Kirshncr 
& Khan 1974, Eastman & Kirshner 1989, Eastman et 
al. 1996, Schmidt et al. 1994ab, Leonard et al. 2002ab, 
Hamuy & Pinto 2002) and the Spectral-fitting Expanding 
Atmosphere Method (SEAM; Baron et al. 1995, 2004), 
to provide distance estimates for objects even within the 
Hubble flow. However, such distance determinations will 
benefit from a thorough modeling of the SN spectrum and 
light curve, as well as from an understanding of the physics 
of spectrum formation. 

Model atmospheres and techniques, which have been 
tested so heavily and applied so successfully to hot 
luminous stars, constitute excellent tools for such SN 
analyses. This is being done with a number of codes: 
WM-BASIC (Pauldrach et al. 2001, Pauldrach et al. 
1996) PHOENIX (Hauschildt 1992ab; Hauschildt & Baron 
1995), EDDINGTON (Eastman & Pinto 1993). Our anal- 
yses are based on CMFGEN (Hillier & Miller 1998), 
which has been used to model OB stars, Luminous Blue 
Variables (LBVs) and Wolf-Rayet (WR) stars (Hillier et 
al. 1998, Dessart et al. 2000, Crowther et al. 2002, Hillier 
et al. 2003). 

O star winds are weak and scattering dominated, quite 
like SN outflows. In the weakest- wind O-stars, line forma- 



tion occurs over a small spatial region so that sphericity 
effects are unimportant. On the contrary O-supergiants, 
where sphericity effects are important, show a range of 
line formation processes, with both photospheric and wind 
sites of emission, absorption and scattering (Hillier et 
al. 2003, Martins et al. 2004, Crowther et al. 2002). 
Conversely, LBV winds are very extended, dense and rel- 
atively cool, and have an ionization structure that resem- 
bles quite closely that of SN. Processes such as Rayleigh 
scattering, charge exchange and two-photon decay are of 
relevance in both LBV and SN outflows (Hillier et al. 1998, 
Hillier et al. 2001, Najarro et al. 1997). WR winds are 
dense and fast, and typically show a much higher ioniza- 
tion than SN. Their composition is typical of progenitors 
of Type II (H-rich) and Type Ib/c (H-free) SN (Dessart 
et al. 2000, Crowther et al. 2002). 

Despite the similarities, SN conditions depart from 
those of hot star winds in a number of ways. The flow is ac- 
celerated only temporarily, say for about a day after core- 
collapse or explosion, followed by a homologous expansion 
in which the velocity follows a Hubble law [v/r constant 
- Woosley 1988). Consequently the outflow density profile 
does not follow an inverse square law as imposed by mass 
conservation (Castor, Abbott & Klein 1975), but instead 
results from the disruption of the envelope by the shock 
wave arising from the explosion (Woosley 1988). Indeed 
models of such explosions, as well as spectroscopic anal- 
yses, support a very fast decreasing density. This rapid 
density decline, together with the large velocity gradient, 
are key ingredients controlling the escape of radiation from 
the flow. 

Another key component for the modeling of outflows is 
the origin of opacity. In H-rich environments, such as mas- 
sive star winds and Type II SN, the dominant source of 
continuum opacity is electron scattering. Although elec- 
tron scattering is more important in Type II SN, there 
is still significant true absorption, due to bound-free and 
free-free processes, to allow the diffusion approximation, 
and the assumption of LTE, to be applied at the base of 
the outflow. For H-free conditions, bound-free and bound- 
bound transitions of metals will become the dominant 
source of opacity. There, including as many levels and 
lines is mandatory for recovering LTE and a thermalisa- 
tion layer at the base. 

Because of their closer similarity to hot massive stars, 
we concentrate on Type II SN in this first series of papers. 
In this first paper, we present the necessary adjustments 
made to CMFGEN in order to model the emergent spec- 
tra from SN outflows, and provide some key insights into 
SN spectrum formation. In Sect. [2 we present the code of 
Hillier & Miller (1998), giving the necessary background 
and highlighting key differences between CMFGEN and 
other models developed and used by different groups. We 
also describe the adjustments from hot star to SN condi- 
tions. In Sect. |31 we present typical results of our spectro- 
scopic modeling, selecting a sequence of observations taken 
during the photospheric-phase evolution of SN1987A and 
SN1999em. These serve as illustrations for the model be- 
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havior, as well as a basis for our interpretation of key 
features of Type II SN during the photospheric phase. We 
also give model parameters, discuss outflow properties and 
provide line identifications. In Sect.^J based on the mod- 
els shown in Sect. |3 as well as the large parameter space 
investigated for the preparation of this paper, we describe 
generic features of Type II SN and how these are sen- 
sitive to model parameters. We discuss the handling of 
the H-recombination front in cool models (Sect. EflT|l . the 
potential problem with the choice of outer grid location 
(Sect. the importance of line-blanketing fSect. 1131) . 
the influence of the luminosity (Sect. ^3} an d density ex- 
ponent (Sect. 14.51) on the spectral morphology, and the 
modeling of hydrogen, helium and nitrogen line profiles 
in the optical range (Sect. H2J|. In Sect. [31 we describe in 
detail the process of line formation in Type II SN, giv- 
ing in particular a physical explanation for the observed 
blue-shift of Balmer lines in photospheric-phase spectra 
of Type II SN. In Sect.® we present our conclusions and 
map out our future endeavors. 

2. Presentation of the code 

Our effort to analyze the light emanating from SN outflows 
makes use of the model atmosphere code CMFGEN, dis- 
cussed in detail by Hillier & Miller (1998, 1999). CMFGEN 
solves the radiative transfer equation in the comoving 
frame, subject to the constraints of radiative and statisti- 
cal equilibrium. The outflow is assumed to be steady-state 
and to have a spherical symmetry. The general procedure 
for analyzing a given object is to input global character- 
istic quantities such as luminosity, radius, chemical com- 
position, velocity and density. 

In hot star outflows, the velocity law is taken from de- 
tailed hydrodynamical computations of radiatively-driven 
winds. For SN outflows, after about one day after the 
explosion, the outflow no longer accelerates but instead 
freely-expands into the surrounding medium. Such an ho- 
mologous expansion is characterized by a Hubble velocity 
law, v(r) = vq (r/Ro) where Rq is the base radius where 
the luminosity is input, and vq the velocity of this layer. 
This layer in Type II SN is chosen such that it is opti- 
cally thick at all wavelengths and thus corresponds to a 
thcrmalisation layer where photons have a Planckian dis- 
tribution at the local temperature T = (i*/47rcri?§) 1 / 4 , 
a being Stefan's constant. Note that this lies significantly 
below the photosphere which we define as the layer where 
the total continuum optical-depth is 2/3 (this includes 
bound-free and free-free processes from all species plus 
electron-scattering), located typically around 1.3-1.5 i?o- 
As time proceeds, (constant-velocity) mass shells expand 
and become more optically thin, so that in the frame co- 
moving with a given mass shell, the thermalisation layer 
recedes deeper into the flow. At present, the time evolu- 
tion is taken into account by fitting not only the spec- 
tral energy distribution (SED) but also the absolute flux 
level at all times. The former gives a tight constraint on 
the ionization balance of the flow and the temperature of 



the thermalisation layer. The second constrains the lumi- 
nosity (or radius of the thermalisation layer). In practice, 
we vary the outflow ionization and the absolute flux level 
by adjusting the radius scale of the model, controlled by 
the parameter Rq, and the luminosity (or the flux) go- 
ing through this inner model boundary, controlled by the 
parameter L». 

For the distance, which is required for the luminosity 
constraint, we have adopted the Cepheid or EPM distance 
to the object under investigation. As discussed later in 
this paper, for a given outflow ionization, the SED is only 
sensitive to large changes in luminosity so uncertainties in 
the distance, at the 50% level, will have negligible effects 
on the SED. 

We adopt an analytical density law and adjust it at 
different epochs according to whatever nature demands. 
Later we show that the steepness of the density fall- 
off is mostly constrained by P-Cygni line profile shapes. 
The analytical form adopted is a power law, i.e. p(r) = 
pa(Ro/r) n . Former investigations, both of the explosion 
itself as well as others on the SED have shown that n is 
large, with a nominal value of 10, and tends to decrease 
over time (see, e.g., Arnett 1988, Lucy 1987, Eastman & 
Kirshncr 1989). po is adjusted so that the Rosseland-mean 
or electron-scattering optical depth at the base is of the 
order of 50. This ensures photons are thcrmalised in the 
base layer. 

In radiatively-driven winds of hot stars, the density 
is usually taken to be clumped on small scales since this 
is expected both theoretically (Owocki & Rybicki 1984, 
Dessart & Owocki 2003) and observationally (Robert 
1992). In multi-dimensional simulations of SN explosions 
(Kifonidis et al. 2003, Reinecke et al. 2002), density con- 
trasts in the heterogeneous outflows are of the order of a 
few at most, thus more than an order of magnitude less 
than suspected in hot stars. Hence, we assume no clump- 
ing of the outflow in our simulations. As we show in the 
next section, the quality of the fits is excellent, suggesting 
that clumping, if present, remains small. 

In addition to clumping, SN outflows are expected 
to exhibit turbulence on a variety of scales (see, e.g., 
Kifonidis et al. 2003, Burrows et al. 1995). We generally 
adopt a microturbulent velocity of 100 km s -1 , which is 
just a few percents of the photospheric velocity - this has 
the effect of broadening the intrinsic (Gaussian) line pro- 
file width, increasing the probability of line-overlap but re- 
ducing the required frequency sampling for the radiative- 
transfer solution. Large scale turbulent motions are ig- 
nored. 

In Type II SN, only a small amount of metals is nucleo- 
synthesized during the explosion (Woosley & Weaver 
1995), so that it seems suitable to assume strict radiative- 
equilibrium. To ensure that such a situation prevails, we 
limit our investigation to early times during the photo- 
spheric phase, say the first month after maximum B band 
light. 

A key asset of CMFGEN is that all lines of all 
species are treated in non-LTE. Since adequate model- 
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ing of the blanketing effects of iron are essential, and 
since such a metal has a very complex atomic struc- 
tures with lots of lines, the atomic structure input in the 
code is re-organized in terms of "super-levels" (SL). In 
this scheme nearby levels, having similar departure coeffi- 
cients, are grouped together and treated as an individual 
level. Despite this grouping, all lines are treated at their 
correct wavelengths (ignoring possible difficulties with the 
atomic data). Applying this procedure to all species per- 
mits the treatment of many more ions and bound-bound 
transitions. In our approach we do not need to assign an 
absorption/scattering cross section to the line blanketing. 
This is automatically determined by the coupled solution 
of the transfer and statistical equilibrium processes. If the 
SL assignment is too coarse, branching (i.e., the decay of 
the upper level to alternate states) may not be treated 
correctly. In particular, to treat the Bowen fluorescence 
mechanism, special care must be taken when making the 
SL assignments. However, numerous tests, in a wide vari- 
ety of conditions, have shown that the SL approach can 
yield an accurate temperature and ionization structure of 
the envelope. Importantly, it is relatively easy to change 
the SL assignments to test their effect. 

We show in Table 1 the largest atomic model used 
in our computations, showing the number of "full-" and 
"super-" levels, and the corresponding number of transi- 
tions. For coarse analyses, only H, He, C, N, O and Fe are 
included. In fact, including Fell ensures that the overall 
temperature structure and ionizing flux are not signifi- 
cantly affected when performing the more detailed com- 
putation with the larger model atom. The atomic data 
come from a wide variety of sources, the Opacity Project 
(Seaton 1987; The Opacity Project Team 1995, 1997), the 
Iron Project (Pradhan ct al. 1996; Hummer et al. 1993), 
Kurucz (1995), and the Atomic Spectra Database at NIST 
Physical Laboratory being the principal sources. Much of 
the Kurucz data was obtained directly from the Center 
for Astrophsyics (Kurucz 1988, 2002). Individual sources 
of atomic data include the following: Bautista & Pradhan 
(1997), Becker & Butler (1995), Butler et al. (1993), Fuhr, 
Martin, & Wiese (1988), Kingdon & Ferland (1996), Luo 
& Pradhan (1989), Luo et al. (1989), Mendoza (1983), 
Mendoza et al. (1995), Nahar (1995, 1996), Nahar & 
Pradhan (1996), Neufeld & Dalgarno (1987), Nussbaumcr 
& Storey (1983, 1984), Peach, Saraph, & Seaton (1988), 
Storey (1988), Tully, Seaton & Bcrrington (1990), Wiese, 
Smith, & Glcnnon (1966), Wiese, Smith, & Miles (1969), 
Zhang & Pradhan (1995, 1997). Note that at present, we 
are missing some species for the modeling of late-stage 
spectra, most importantly Ci, Sen and Ban, something 
we will remedy in the future. 

Before describing results obtained with CMFGEN, it 
is worthwhile stressing the differences with other models 
used for SN studies. As described in the introduction, a 
number of general purpose model atmosphere codes, often 
developed originally for the treatment of the expanding 
outflows of hot stars, have been applied to SN. The general 
practice for the study of Type II SN is to assume spherical 



Table 1. Summary of the model atom used in our calcu- 
lations. Np is the number of full levels, iVg the number of 
super levels and At the corresponding number of transi- 
tions. 
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symmetry and steady-state, although some attempts are 
now made to perform 3D Monte-Carlo radiative-transfer 
calculations (Thomas et al. 2003). The level of refinement 
for the radiative transfer computation varies significantly: 
many studies treat all metal lines in LTE, using non- 
LTE for H and He only (Schmutz et al. 1990; Eastman 
et al. 1989, 1996). Although much more computationally 
expensive, a higher level of consistency is met by stud- 
ies that treat instead all lines in non-LTE (Hoflich 1988; 
Mitchell et al. 2002). Finally, the alternative parameter- 
ized approach with SYNOW, tagged "direct-analysis" , is 
less physically consistent than previous methods, but con- 
stitutes a useful tool to identify spectral features - a chal- 
lenging task in SN spectroscopy - and can thus serve as a 
basis for more accurate model atmosphere computations. 

Before proceeding further it is worth giving an hon- 
est appraisal of some of the deficiencies of the present 
study. First we assume spherical geometry and homogene- 
ity. Many spectropolarimetric studies have shown that at 
least for some SN, a departure from spherical symmetry on 
a large scale is present, although mostly restricted to late 
times (see, e.g., Kasen et al. 2003 for Type la SN, Hoflich 
et al. 1996 for Type lb SN, and Leonard et al. 2001 for 
Type II SN). Further, due to Rayleigh- Taylor instabilities, 
there is likely to be substantial mixing of material from the 
inner layers into the outer layers (Kifonidis et al. 2003), 
as well as smaller scale structures (with a scale of the or- 
der of or below the photospheric radius value). One such 
evidence for inhomogencities was seen in Ha line profile 
variations in late-time spectra of SN 1987A (Spyromilio et 
al. 2003). Second, we neglect the time dependence of the 
flow. This is certainly not valid for modeling the luminos- 
ity evolution from hydrodynamical calculations, especially 
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at early times, but is expected to be a reasonable approxi- 
mation in the atmosphere where we can fix the luminosity 
using observations. 

At present we ignore relativistic effects other than the 
frequency shift arising from the first order Doppler cor- 
rection. Such effects have been discussed by Hauschildt et 
al. (1991), and will be the subject of another paper. All 
these effects are crucial for hydrodynamical modeling, but 
less so in the approach we adopt where we constrain the 
emitted luminosity by the observations, and only treat the 
"photosphcric" layers. Other effects, such as an accurate 
treatment of line blanketing, are much more important in 
determining the observed spectrum. The primary motiva- 
tion for doing so was to utilize CMFGEN in a standard 
configuration — this allowed convergence difficulties to 
be attributable to the properties of SN envelopes, rather 
than complications arising from the relativistic transfer. 
In an earlier attempt at modeling SN the strongly scat- 
tering envelope gave rise to difficulties in the temperature 
convergence which were wrongly attributed to problems 
in the relativistic transfer, whereas they were related to 
the use of SL (Hillier 2003). 

Due to the large photospheric velocities, we modified 
the gray temperature formulation to take into account the 
first order Doppler shifts. This formulation was necessary 
to provide an improved estimate of the temperature struc- 
ture at depth at the beginning of the model iteration. 

3. Spectroscopic modeling 

In the present paper, we wish to discuss the basic proper- 
ties of a few models whose synthetic spectrum fits reason- 
ably well the observations at early, intermediate and late 
times of the photosphcric-phasc evolution of Type II SN. 
We base our discussion on two very well observed Type II 
SN, SN1987A and SN1999cm. The former was associated 
with the neutrino detection on February 23.316 UT 1987 
with the Kamiokande detector (Koshiba 1987). Here, we 
use the first spectrum taken both in the UV and in the 
optical (Pun et al. 1995, Phillips et al. 1988), on February 
24th, thus only one day after discovery (and collapse). We 
use this observation for the early-stage description of Type 
II SN spectra. For SN1999em, no neutrino detection was 
made and thus no precise dating of the core-collapse of 
the progenitor exists. Inferences, based on the light curve 
appearance or the EPM, suggest the explosion took place 
at around October 27th 1999, with an uncertainty of 1-2 
days. It was absent in images of the host galaxy NGC 1637 
taken on October 20.45 UT, while first detection of the SN 
was made on October 29th (JD 2451480.94; Li 1999). The 
resemblance of the spectrum taken on October 30th with 
the spectrum of SN1987A taken on the 24th of February 
1987, indicate that SN1999em was indeed caught at a 
very early phase of evolution. Here, we use observations 
of SN1999cm taken on the 5th of November 1999 (ca. one 
week after explosion; Baron et al. 2000, Leonard et al. 
2002a) and on the 14th of November 1999 (ca. two weeks 
after explosion: Hamuy et al. 2001) to describe respec- 



tively the typical intermediate- and late-stage appearance 
of Type II SN spectra. Such time qualifications to obser- 
vations are to be taken in a broad context. We defer until 
a follow-up paper a detailed analysis with a fine time- 
sampling of the evolution of the SED of a Type II SN 
during the entire photospheric phase. 

Let us describe the spectral evolution in the 
UV/optical ranges of these two well observed SN. For 
each SN, the first spectrum taken shares striking similar- 
ities: the flux in the UV is much larger than in the opti- 
cal where one sees hydrogen Balmer lines and He 1 5875A. 
Following this stage, as the outflow expands and cools, 
the UV flux diminishes and the He 1 5875 A P-Cygni pro- 
file weakens. Depending on the object, this takes a range 
of time scales: SN1987A did that step in just one day, 
while the corresponding spectral evolution took 10 days 
for SN1999em. One can then identify a cooler stage where 
both the UV flux and He I lines are completely gone. Metal 
lines start appearing in the optical range and the Can 
multiplct around 8500Ais present. Further, the Nai dou- 
blet at 5890A starts appearing as a weak P-Cygni pro- 
file. During this entire sequence of events, hydrogen re- 
combination becomes more and more pronounced so that 
the next stage is characterized by the presence of the hy- 
drogen recombination front close to the photosphere. The 
SED (i.e., f\) then peaks at or beyond 4000A. Subsequent 
evolution is more gradual, with significant changes taking 
place over longer time-scales. 

When fitting synthetic spectra to observations, we red- 
den the theoretical SED using the Cardelli law (Cardelli 
et al. 1988) for SN1999em and the Seaton (1979) and 
Howarth (1983) laws for the Galaxy and the LMC re- 
spectively. For SN1987A, we take the Cepheid distance 
to the LMC of 52kpc for SN1987A and adopt a redden- 
ing value E(B - V) = 0.05 for both galactic and LMC 
components. This is slightly less than the value 0.16 pro- 
posed by Lundqvist & Fransson (1996) but seems war- 
ranted to reproduce the observations of SN1987A on the 
24th of February. For later times, this choice has essen- 
tially no effect on the fit quality to observations. For 
SN1999em, we use the Cepheid-based distance of 11.7Mpc 
to NGC 1637, somewhat different from the EPM based dis- 
tance of 8.2Mpc (Leonard et al. 2002a) or 7.5Mpc (Hamuy 
et al. 2001). These differ slightly but within such limits, 
the computed SED remains identical, only the absolute 
flux level to reproduce observations will change, by about 
a factor of two. Note finally that we correct for the re- 
cession velocity of the galaxy host of SN1999em by blue- 
shifting the observations by 770 kms -1 (Leonard et al. 
2002a, Baron et al. 2000). 

Due to the ubiquitous presence of line-overlap, it is in- 
structive to investigate which lines contribute to observed 
"line features". Except on very rare occasions and/or 
for very few lines (early-spectra in the optical, Ha for 
most epochs), the very large expansion velocities are such 
that observed features are composed of different lines of 
different species. If these arise from species with differ- 
ing ionization potential, they are likely to form in dis- 
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SN1987a, 24/02/1987, E(B-V) = 0.05 (Gal), E(B-V) = 0.05 (LMC) 





I i I i i I I i i i I i i i I i 

2000 4000 6000 8000 

A (A) 

Fig. 1. Top: Illustration showing the comparison between reddened model spectra (blue) and observations of SN1987A 
(black) taken on the 24th of February 1987 (Pun et al. 1995, Phillips et al. 1988). Each model flux distribution is 
scaled for an assumed LMC distance of 52kpc. In red we show the theoretical continuum energy distribution of the 
model. Model characteristics are: L* = 3 x 10 8 L©, i? p hot = 2.74 x 10 14 cm (or 3940 R Q ), v phot = 17, 700kms -1 , 
T r ff = 11.200K, n = 12, p n hot = 1-1 x 10 -13 gem -3 . Bottom: Montage of rectified spectra computed by including 
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tinct regions of the outflow and therefore have a profile 
width and strength quite unpredictable a-priori. To assess 
the importance of this issue, our code can compute the 
emergent spectrum from a formal solution of the trans- 
fer equation based on the full non-LTE solution, includ- 
ing all continuum processes (all bound-free and free-free 
processes) but including only the line transitions of the 
desired species/ionization stage. Although this then ig- 
nores the non-linear interaction of overlapping lines on 
the observed spectrum, it allows us to assess with enough 
accuracy what species contribute to a given feature. 

Finally, in each of the following sections, we give details 
on model parameters corresponding to the synthetic spec- 
tra shown in the figures. However, unless otherwise stated, 
all simulations use a chemical composition adequate for 
a young supergiant massive star. In practice, we use the 
following CNO-cycle equilibrium values for a massive pro- 
genitor star (Prantzos et al. 1986): slightly enriched abun- 
dances of helium and nitrogen H/He = 5, N/He = 6.8 x 
1CT 3 , and depletion of C and O, i.e. C/He = 1.7 x 1(T 4 
and O/He = 1CP 4 (all given by number). Additionally, 
when the model refers to SN1999em, the abundance of 
metal species (excluding CNO elements) is set to its Solar 
Neighborhood. This is compatible with the inferred so- 
lar or over-solar metallicity of its host galaxy NGC 1637 
(Leonard et al. 2002a). For SN1987A, metal abundances 
are all scaled by a fixed factor of 0.4, compared to their 
solar value. 



3.1. Early stage 



Early— stage mode 




Fig. 2. Relative ionization for hydrogen (solid), helium 
(dashed) and iron (dash-dotted line) for the early-stage 
model of Sect. 13. ll and Fig.^ The ionization fractions are 
expressed relative to the abundance of all species. Also 
shown is the electron density on the top axis [color]. 

We show in the top panel of Fig. synthetic fit 
(color) to observations (black) of SN1987A taken on 
the 24th of February 1987. The red curve represents 
the continuum energy distribution of the model, while 



the blue curve corresponds to the synthetic flux distri- 
bution when all line and continuum processes are in- 
cluded. The model characteristics are the following: we 
have L» = 3 x 10 8 Lo, i? phot = 2.74 x 10 14 cm (or 
3940i? Q ), v phot = 17, 700 kms- 1 , T cff = 11,200K, n = 12, 
Pphot = 1-1 x 10 -13 gem -3 . This model is computed using 
the full model atom described in Table 1. We had trouble 
fitting both UV and optical spectral regions, with the ad- 
ditional constraint posed by line profiles. Here, (enhanced) 
UV flux was partly achieved by adopting a slightly lower 
reddening than usually employed (0.1 compared to 0.16). 
Obtaining weak absorption and emission in P-Cygni pro- 
files could only be done by choosing a high density expo- 
nent, say above 8, while much higher values (of the order 
of 20) led to very weak and anomalous P-Cygni absorp- 
tion. Density exponent values of 10-12 seem best suited 
for optical fits to line profiles in this spectrum, somewhat 
higher than the value of 6-7 proposed by Lucy (1987). 
Admittedly, our fit to the Ha trough and H/3 peak flux 
could be improved (see inserted box in top-right corner 
in Fig. nj ; a more gradual density fall-off in the regions 
where these lines form would likely strengthen the hy- 
drogen Balmer lines and extend the associated P-Cygni 
troughs (see Sect. JSJ). 

The outflow ionization is high with a SED that peaks 
in the UV. The line-blanketing due to metal species is 
modest — a significant UV flux is still present. We show 
in the bottom panel of Fig.^the actual contribution from 
each of the species present, by computing the synthetic 
spectrum including only the bound-bound transitions of 
a given species. The resulting SED is then normalized by 
the continuum energy distribution and ordered in Fig. ^ 
by increasing atomic weight, starting from the bottom. 
To enhance the visibility of weak optical features, that 
nonetheless contribute to observed features, we apply a 
scaling to the rectified spectra so that we in fact show 
(fx/ fc — 1) x xscale + 1, with xscale = 4 beyond 3500A, 
and unity elsewhere. 

Few features are actually associated with a small 
set of lines. At 2800A, the feature is due to Mgn(3p- 
2s) 2802. 7A and the doublet Mgn(4s-3p) 2928.6-2936.5A. 
At 1900A, the feature is due to a small but larger num- 
ber of Aim lines, predominantly Alni(3p-3s) 1859A and 
(4f-3d) 1935A. Other UV features are the result of a large 
number of lines, in particular from Fe in and Ni III. Indeed, 
the combination of a rich line-spectrum of metals, partic- 
ularly that of iron, together with the large Dopplcr shifts 
induce a strong line-blocking effect on the radiation field. 

In the optical (see the inserted zoom in the top right 
corner of Fig.^l, we find the usual suspects of Type II SN, 
i.e. hydrogen Balmer lines: Ha at 6562. 8A, H/3 at 4861. 3A, 
H7 at 4340. 5A, RS at 4101.7A, and He at 3970.0 A. At 
shorter wavelengths, the higher transitions in the Balmer 
series merge together until the Balmer jump at 3646A. 
The later region might be weakly contaminated by the 
presence of the Ca 11 doublet (4p-4s) at 3933.6-3968.5A. 

We can also find the presence of the Hei5875A. For 
SN1987A, this line is only present in the first spectrum, 



8 Luc Dessart, and D. John Hillier: Type II supernovae spectroscopic modeling 



SN1999em - 5th of Nov. 1999 - E(B-V) = 0.1 
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Fig. 3. Top panel: Synthetic fit (color) to HST and Lick observations (black) of SN1999em taken on the 5th 
of November 1999 (Baron et al. 2000, Leonard et al. 2002a). The red (blue) curve corresponds to a model with 
L* = 4xl0 8 £ Q (L* = 5xl0 8 L Q ) andTcff = 8,850K (T cff = 9,200K). Other properties are identical: i? ph ot = 6.6 xlO 14 
cm (or 9540i?©), w p hot = 8, 750kms , n = 10 and /9 p hot = 4.1 x 10~ 14 gem -3 . Lower panel: Montage of rectified 
spectra (same as in Fig. IB [color 1 . 



Luc Dessart, and D. John Hillier: Type II supernovae spectroscopic modeling 



9 



shown here, while for SN1999em, this line is observed 
for about a week (see Sect. 14.6.21) . It serves as a key di- 
agnostic of the outflow ionization since for just slightly 
cooler temperatures, it becomes optically-thin and its 
strength vanishingly small. Hei5875A is the one helium 
line that gets most attention because it is isolated and 
thus, uncontaminated. However, the model predicts the 
presence of a group of Hei lines around 4000 A (Hel(3p- 
2s) 3888.6A, He i(4p-2s) 3964.7A, He i(7d-2p) 4009.2A, 
Hei(7s-2p)4024.oA, Hei(5d-2p)4026.2A) as well as 
around 4500A (Hei(5d-2p) 4388-lA, He i(5s-2p) 4437.5A, 
Hel(4d-2p)4471.5A, He i(4s-2p) 4713. lA), followed by 
Hei(4d-2p)4921.9A and Hei(3d-2p) 6678A. The strength 
of the resulting features is predicted to be of the same or- 
der as that of He 1 5875A and hence, it is clear that helium 
lines are present throughout the optical range. 

We also find the presence of Nil lines in the opti- 
cal: in the blue wing of He 1 5875A as well as in the 
blue wing of H/3, at respective wavelengths of 5495-5680A 
and 4477-4630A (each of these is part of a multiplet). In 
SN1987A, these features are quite weak (resulting to a 
large extent from its unusually high expansion velocity) 
but their unambiguous presence is seen in the early spec- 
tra of SN1999em (see Sect. I4.fi.3jl and even more so in 
SN1999gi (Leonard et al. 2002b). We defer until Sect. |4~Q 
for a more detailed discussion on these N II lines and their 
presence in early spectra of Type II SN. For carbon and 
oxygen, their expected under-abundance during the CNO 
cycle due to mass loss (Prantzos et al. 1986), which we 
adopt here, is such that not a single line from C II, C III, 
O I, O II, O in is predicted to be noticeable in this model. 
Note also that one can see few weak features arising from 
Si in and Fe in lines in the optical, while most metal fea- 
tures are in the UV at this epoch. 

Finally, we display the outflow ionization in Fig. [5] 
Hydrogen is fully ionized, helium doubly ionized only at 
the base, singly ionized at intermediate heights and neu- 
tral in the outer part of the outflow. Iron is present in 
the form of Fe 4+ inside through to Fe 2+ outside, with Fe + 
never setting in. Note that in this model, the photosphere 
is located at r = 1.48i?o or log(r/i? ) = 0.17, thus in a 
region where H + , He + and Fe 2+ dominate. 

3.2. Intermediate stage 

In this and the following section, we turn to the subse- 
quent evolution of Type II SN spectra, using observations 
performed at ca. 1 and 2 weeks after the estimated ex- 
plosion time of SN1999em (Li 1999; see introduction to 
Sect.©. 

We show in Fig. synthetic fits (color) to observa- 
tions (black) taken on the 5th of November by Baron 
et al. (2000) using the Hubble Space Telescope, as well 
as optical observations from Leonard et al. (2002a). The 
two models differ only in luminosity: I, = 4 x 10*^0 for 
the red model and L, = 5 x 1O 8 L0 for the blue model. 
The resulting model properties differ only in effective tem- 
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Fig. 4. Ionization structure for the intermediate-stage 
model described in Sect. 13.21 and shown in Fig. EI The ion- 
ization fractions are expressed relative to the abundance 
of all species [color]. 



perature: T eff = 8, 850K (red model) and T off = 9, 200K 
(blue model) while the common properties (parameters) 
are the following: i? p hot = 6.6 x 10 14 cm (or 9540i?©), 
Wphot = 8,750kms _1 , n = 10 and p phot = 4.1 x 10~ 14 
gcm~ 3 . The red model would give a better fit to obser- 
vations if we chose E(B — V) = 0.08 — 0.09 instead of 
the adopted value of 0.1. This shows that the UV range is 
important to place very strong constraints on the model 
temperature as well as the adopted reddening, helping to 
lift the degeneracy in the impact of these two parameters. 

In the bottom panel of Fig. |3| we show the synthetic 
spectrum for the blue model when only the bound-bound 
transitions of selected species (indicated on the right) are 
included together with all bound-free and free-free pro- 
cesses. As in Fig. ^ we apply a scaling of a factor of four 
to the rectified spectra beyond 3500A to enhance the vis- 
ibility of weak features. 

At this stage of its evolution, the SN outflow has cooled 
significantly, so that Fe n line-blanketing starts to strongly 
diminish the UV flux. This blanketing is not so severe, so 
that some residual UV flux is observable, carrying infor- 
mation on the species at its origin. From the lower panel, 
we see that the line-blanketing makes the UV flux smaller 
than the continuum flux when all lines are ignored. In 
other words, emission features in the UV are to be un- 
derstood as spectral windows where this line-blanketing 
is less relative to nearby regions, rather than an excess 
flux above the continuum level. We thus identify enhanced 
line-blocking due to Fen, Fern, Nin, Sin, Alii, Aim. An 
exception to this is due to Mgn lines, already discussed 
in the previous section, which show pronounced emission 
and absorption around 2800A. 

In the optical range, besides hydrogen lines, Hei fea- 
tures are numerous with a strength of about one fifth 
of Ha (see previous section for a census of these lines). 
However, the outflow has cooled down too significantly to 
show N II lines, and we instead see the appearance of weak 
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SN1999em - 14/11/1999 - E(B-V) = 0.1 
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Fig. 5. Top panel: Synthetic fit (blue) to observations (black) of SN1999em taken on the 14th of November 1999 
(Leonard et al. 2002a) . We also show the continuum flux level (red) computed by ignoring all bound-bound transitions in 
the formal solution of the radiative transfer equation. The model parameters are: L» = 1.5 x 10 8 L©, i? p hot = 6.15 x 10 14 
cm (or 8840i?o), u p hot = 6, 350 kms -1 , T e g = 6, 800K, n — 10, p p hot = 8.7 x 10~ 14 gcm~ 3 . In the top panel, we insert 
a zoom on the Nai5596-5590A region: the green curve corresponds to the same model as above (blue) but with the 
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metal lines in the optical spectrum. We have the Mgn(4f- 
3d) triplet at 4481 A. We also have Sin(5p-3d)4076.8, 
Siii(4f-3d)4128.0A, the Sin(4d-4p) doublet at 5056.0A 
and the Sill(4p-4s) doublet at 6347.lA and 6371.4A. We 
finally attribute some features to Fell and Fein but due 
to the large number of contributors, we cannot enumerate 
them. Note nonetheless that contrary to identifications by 
Leonard et al. (2002a), the feature around 5200A has Fe III 
lines (groups of lines around 5080A and 5 160 A) as main 
contributors, rather than Fe II, and also overlaps with lines 
of Sin at 5056A and Hei at 4921.9A and H/3. This careful 
log shows that it is difficult to associate "features" in SN 
spectra with individual lines of selected species. 

We show in Fig.0|the outflow ionization for hydrogen, 
helium and iron for the blue-curve model of Fig. It re- 
sembles quite closely the previous case. However, the He + 
region has shrunk, causing Hei lines to weaken. Note also 
that H 0+ and Fe + have a larger fractional abundance in 
the outer part. The photosphere is located at r — 1.59i?o 
(log(r/i?o) = 0.2), where H + and Fe 2+ dominate while 
helium is roughly equally present in its first ionized and 
neutral forms. 

3.3. Late stage 

We now turn to the description of a Type II SN spectrum, 
corresponding to a late-stage in the photosphcric-phase 
evolution. In Fig. [5] we show synthetic fits (color) to the 
observations (black) of SN1999em taken on the 14th of 
November 1999. The red curve corresponds to the con- 
tinuum energy distribution of the model, while the full 
synthetic spectrum is shown in blue. The model parame- 
ters are: L* = 1.5 x 10 8 L Q , i? p hot = 6.15 x 10 14 cm (or 
8840i? Q ), v phot = 6,350kms-\ T cff = 6,800K, n = 10, 
Pphot = 8.7 x 10~ 14 gem" 3 . 

Only optical data is available for SN1999em at this 
date so we rely entirely on the optical range to assess the 
UV flux level. The absence of He I lines and the increasing 
strength of Fe II lines are consistent with a much cooler ef- 
fective temperature compared to the previously discussed 
models, associated with a strong line-blanketing in the UV 
range. The lower panel shows that this line-blocking stems 
mostly from Fen, and to a lesser extent from Nin, Tin 
and Mgn (the 2800 A feature). We find that as long as a 
big model atom for Fe II is included in our computations, 
the gross properties of the outflow and the emergent spec- 
trum remain unchanged when adding extra metal species, 
e.g. Nin. 

The optical range clearly shows the presence of P- 
Cygni profiles from Can 3933-3968A and Can 8498- 
8542-8662A, which gives further support for a relatively 
low ionization of the outflow. Blends of H7, Tin and Fen 
contribute to the feature at 4300A while blends of Till, 
Mgn and Fen contribute to the feature at 4600A. Around 
5000A we see a myriad of Fen lines which overlap with 
H/3, while Fein lines are no longer predicted. The flux 
goes below the continuum level around 5500A, which we 



interpret as the appearance of the doublet line Na 1 5896- 
5890A. Indeed, taking the same model and enhancing the 
sodium abundance by a factor of four (Fig. top-panel, 
green curve in the inserted box), as expected for a su- 
pergiant progenitor (Prantzos et al. 1986), gives a much 
better fit to the observed sodium doublet. The dip on the 
blue-side of this sodium doublet has been associated with 
a number of Sen lines (Leonard et al. 2002a), a species 
not yet included in our atomic database. Further to the 
red we reach a spectral region affected by weak Fe 11 lines 
which, in combination with Sin6347-637lA gives rise to 
some filling-in of the P-Cygni trough of Ha. 

As before, we show the outflow ionization predicted by 
the model (Fig. - note that the highest ionization state 
of iron shown is Fe 3+ but the lowest one is now Fe 0+ ). 
Unlike the previous two cases hydrogen recombines to its 
neutral state, while in the corresponding region, iron is 
present in its Fe + state (giving rise to the Fe II lines) . Note 
that in the regions where hydrogen recombination takes 
place, the electron-density goes down more steeply with 
radius (this is seen by inspection of the spacing between 
tick marks on the top axis). The photosphere is located at 
r = 1.48i?o (l°g( r /-Ro) = 0.17), right in the recombination 
front: hydrogen (iron) is equally represented in H + and 
H 0+ forms (Fe 2+ and Fe + ) while helium is fully neutral. 
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Fig. 6. Relative ionization for hydrogen (solid), helium 
(dashed) and iron (dash-dotted) for the late-stage model 
of Sect. 13.31 Note that the highest ionization state of iron 
shown is Fe 3+ rather than Fe 4+ , but the lowest one is now 
Fe 0+ rather than Fe + [color]. 



4. Discussion on model results and properties 

4.1. Hydrogen-recombination front in cool models 

For models having an ionization similar or lower to the 
case shown in Sect. 13.31 (late-stage model), we had to in- 
troduce an adaptive grid in CMFGEN in order to resolve 
the region where hydrogen recombines. Indeed, this tends 
to occur over very restricted spatial scales in the outflow. 
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We show in Fig. [7| a plot of the outflow ionization for 
a very cool model with T c g = 6000K. Note how steep the 
recombination of hydrogen is as well as the concomitant 
recombination of iron from Fe 2+ to Fe + . Also shown is 
the radial variation of the electron-density (top axis); tick 
marks are more closely packed around the front, following 
the strong variation of the mean-electron number through 
this region. 

We find that the shape of the front depends on the spa- 
tial scale of the outflow, i.e. the bigger Ro the more grad- 
ual the front is. Further, for identical model parameters, 
a smaller turbulent velocity gives reduced line-blanketing 
and thus a more extended hydrogen recombination front. 
Reducing the hydrogen abundance also reduces somewhat 
the steepness of the front. This may result from the fact 
that the electron-scattering optical depth is less tied to 
the front itself, which sets the position of the photosphere 
for H-rich models. Finally, varying the density exponent 
(see below) also modifies the shape of the recombination 
front, by extending it over larger spatial scales for smaller 
values of n. 
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Fig. 7. Relative ionization for hydrogen (solid), helium 
(dashed) and iron (dash-dotted) for a very cool model. 
For clarity, we only show the inner region of the outflow. 
The radial variation of the electron-density is also shown 
on the top axis [color]. 



4.2. Choice of outer boundary 

There is a strong computational incentive to limit the ex- 
tent of the simulated grid. Ideally, the outer grid location 
should be where the medium becomes optically thin at all 
wavelengths. Unfortunately, Hi and Hei bound-free edges 
in the UV, as well as a few resonance lines remain optically 
thick even far above the photosphere. This occurs even 
though the density at r is only (Ro/-R m ax)™ its value at Ro, 
n being of the order of 10. The drawback of stretching the 
computation to large radii is that the maximum velocity 
in the outflow can become very large, i.e. up to one third 
of the speed of light for the modeling, e.g., of the spectrum 



of SN1987A on the 24th of February 1987. In expanding 
outflows, a line with frequency vq at a given wind location 
r can in principle interact with photons emitted at r' with 
a frequency between Vq and Vq(1 + w- (v(r) — v(r)')) where 
n is the unit vector linking the two points of interest. The 
frequency sampling around line center must therefore be 
fine enough to account for potential line-transfer effects. 
But the larger the maximum velocity in the computation 
is, the larger is the frequency array and thus the more 
computationally demanding the model is. 

In general, our simulations extend out to lORo- But 
we have investigated whether the choice of a bigger or a 
smaller maximum radius would produce a noticeable effect 
on the SED. Models suitable for the earlier evolution of 
Type II SN during the photosphcric phase (first week) 
tend to require a higher density exponent than at later 
times (see, e.g., Hoflich 1988 for SN1987A; Eastman & 
Kirshner 1989), minimising the effect resulting from the 
modulation of the location of the maximum radius. For 
this check, we therefore have chosen a model with a density 
exponent of 8. This study is done with a small model atom, 
including only the most abundant species Hi, Hei, Cll, 
Ni, Oi, Fen, and Fern. 

We show in Fig. |S| the SED for models with a maxi- 
mum radius of 2, 3, 6, 10, and 15i?o- For this latter case, 
we have modified the Hubble law by forcing it to smoothly 
plateau beyond lORo to reach at 15Ro only 110% of 
the value it has at 10i?o- Only three curves are seen in 
Fig. [SJ Indeed, we find that the SED for the cases where 
R ma i > 6i?o are indistinguishable (they overlap under the 
black curve), a very desirable property which therefore 
justifies a significant reduction of the maximum radius. 
Out of curiosity, we have also reduced the maximum ra- 
dius down to 3 (red curve) and finally 2Rq (blue curve). 
Interestingly, beyond 3000A, the spectrum computed with 
R max = 3Ro, is indistinguishable from that computed 
with R max = 6Rq. This suggests that the continuum and 
line formation beyond 3000A is confined to the photo- 
sphere, located in these models at L8.Ro- However, the UV 
flux lies significantly above the curve corresponding to the 
model with R max = 6Ro, which results from the neglect of 
line-blocking from the outflow layers between 3 and 6Ro- 
Indeed, metal lines in the UV remain optically thick far 
above the photosphere. Going down to R max = 2i?o, we 
see that the UV flux is even more in excess, for the same 
reason given as above, only magnified by the further re- 
duction in maximum radius. But additionally, hydrogen 
lines (most notably Ha, H/3, H7, Pa) are affected, both 
in emission and absorption. Note that the latter is nearly 
absent in Ha, suggesting its P-Cygni trough forms over 
an extended region. In fact, adopting such a small value 
for -R m ax has a strong influence on the non-LTE solution 
of the problem. For example, we now find that the outer 
radius is optically thick throughout the wavelength range, 
so that there is no photosphere in this model! 

Although the tests with i? max of 2 and 3Ro pre- 
sented only pedagogical material, the adequacy of choos- 
ing i? max = 6i?o in our simulations is a very useful result. 
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Fig. 8. Comparison between emergent spectra computed with identical parameters but differing in the adopted 
maximum radius, chosen at 2 (blue), 3 (red) and 6, 10 and 15i?o (overlapping perfectly to form the black curve). To 
facilitate the comparison of the relative synthetic fluxes for the different models, these have all been normalised to 
unity, at 10,000A (left panel) and 20,000A (right panel). Ro is identical in all cases. Note also that the red and black 
curve overlap perfectly everywhere except in the UV. We employ a small model atom for this investigation, i.e. solely 
Hi, He i, Cn, Ni, Oi, Fen, Fein [color]. 



4.3. Influence of metallicity on the emergent spectrum 

As was shown in Sect. line-blanketing becomes very 
effective in blocking the continuum star light as soon as 
the medium has cooled down sufficiently for Fell to be 
present. The effective temperature at which this transition 
occurs is around 8,000K. Starting from the intermediate- 
stage model of Sect. 13.21 we have investigated what the 
impact would be of reducing or enhancing the metallicity, 
i.e scaling all metal abundances by a factor of a half or 
two. We show the results in Fig. ^\ where we reproduce 
the model of Sect. l3.2l in black, the twice-solar (half-solar) 
metallicity model in red (blue). As expected, the higher 
the metallicity, the more significant the line-blocking, and 
this is seen most clearly in the reduction of the emergent 
flux in the UV. Going from the model with Z = 0.5Zq, 
to Z = Zq and Z = 2Zq, the Rosseland-mean optical 
depth at the base of each model is increased by about 
10%, the photosphere moves outwards by less than 1% of 
the radius in the thermalisation layer, with a decrease in 
effective temperature and an increase in the velocity of 
the same order. Hence, despite the fact that each model 
possesses essentially the same outflow properties, the dif- 
ferent adopted metallicity leads to significant differences 
in the UV spectral morphology. The optical spectral mor- 
phology and in particular the slope of the SED is essen- 



tially unchanged, reflecting the similarity of each model 
properties. The only differences are a somewhat weaker 
Hei5875A line at higher metallicity, perhaps resulting 
from the reduced strength of the ionizing flux, and an 
enhanced strength of Fen/Fein lines around 5200 A and 
4500A. 

It thus seems that the UV SED can be used to con- 
strain the environmental metallicity at early times when 
the line-blocking due to the Fell lines is still moderate. 
At later times, e.g., for the late-stage model discussed in 
Sect. 13.31 the UV flux is negligible and thus placing con- 
straints on the metallicity has to be done from the optical - 
there are no strong iron lines at longer wavelengths. Strong 
line-overlap makes this task difficult. Finally, note that the 
probability of seeing metals nucleosynthesized during the 
explosion increases with time so the primordial metallic- 
ity of the progenitor can be more reliably constrained from 
spectra taken at earlier times. 

4.4. Influence of Luminosity 

When modeling SN there is always an uncertainty in its 
distance. In this section we therefore investigate how the 
synthetic flux distribution changes with luminosity. Note 
that in our current approach the luminosity of the model 
is adjusted so that, for an adopted Cepheid distance, the 
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Fig. 9. Comparison of the emergent flux from models with identical parameters apart from a different adopted 
metallicity. The black curve describes a model with solar metal abundance. The red (blue) shows the resulting emergent 
spectra when all metal abundances are scaled by a factor of two (one half). The model spectra beyond Ha are relatively 
insensitive to changes in metallicity — therefore for ease of comparison we have normalised the fluxes to unity at 7000A 
[color] . 



model flux distribution lies within a factor of 2 of that 
observed. In the future, it will be desirable to modify 
CMFGEN in order to use a SN luminosity computed from 
a hydrodynamical simulation of the explosion. 

We show in Fig. ^] a set of three models differing in 
model luminosity (black curve: L„ — 2.5 x 10 8 L©; blue 
curve: = 2.5 x 10 9 L Q ; red curve: T* = 2.5 x 10 7 L Q ). In 
order to maintain a similar total base continuum optical 
depth and effective temperature, we scale up (down) the 
base radius by \/l0 and scale down (up) the base density 
for the model corresponding to the blue (red) curve. To 
position the three models in the same figure, we display 
the flux of the blue (red) model scaled down (up) by a 
factor of 10. With this scaling the model that has a higher 
flux in the UV has a lower flux in the optical and vice- 
versa, a likely result of the modulation of line-blanketing 
due to the variation in spatial scale of the corresponding 
model. Note that we use a small model atom, including 
only the most abundant species H I, He I, C II, N I, O I, Fe II 
and Fein. 

These various model scalings lead to a very similar 
photospheric velocity, i.e. 6,700 km s -1 (red curve) and 
7,200 km s -1 for both black and blue curves, as can be 
seen from the nearly overlapping positions of absorp- 
tion in line features. As expected from Fig. ^| the ef- 
fective model temperature increases from the red (T c ff = 
7,000K), to the black (T off = 7, 900K) and blue model 



(T ff = 8,400K). The reason for this is that the photo- 
spheric density varies a great deal amongst those three 
models, increasing from high to low luminosity models 
(blue: pphot = 0.88x 10 14 gcm~ 3 ; black: p p hot = 2.8x 
10 14 gcm- 3 ; red: p ph ot = 16.3x 10 14 gcm- 3 ). The higher 
the photospheric density the stronger the line-blanketing 
and the more depleted the UV flux appears. 

Thus, when a large flux difference between a model 
and a given observation exists, a satisfactory fit to the 
observed energy distribution will likely become poor once 
the model and observed flux levels are comparable. We 
find that little or no changes are observable in the syn- 
thetic flux distribution when the luminosity is varied by a 
factor of two, so this sets the error one can make on the 
model luminosity when performing fits to observations. 

4.5. Effects of varying the density exponent 

In our approach, we use an analytical description of the 
outflow density, characterized by a base density po,n and 
a density exponent n so that at a given height r, we have 
p{r) = po,n(R*/r) n ■ To ensure that the diffusion approx- 
imation is fulfilled at the base of the envelope, the base 
density is chosen so that the total continuum optical depth 
of the outflow is approximately 50. The density exponent 
is, however, constrained from observations. We have there- 
fore a lot more flexibility than if we based our calculations 
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Fig. 10. Comparison of the emergent flux from models differing in luminosity, enhanced and decreased by a factor of 
ten compared to the model in black. Also, we scale the radius and the base density up or down by vTO, to maintain 
the same photospheric velocity for all (see text for details). The model atom used is small, including only the most 
abundant species, i.e. Hi, Hei, Cli, Ni, Oi, Fen and Fein [color]. 



on hydrodynamical inputs, for which the density distribu- 
tion is set once and for all at the start of the homologous 
expansion. This flexibility is not necessarily a benefit since 
varying the density exponent can influence the emergent 
spectrum in a variety of ways, widening considerably the 
already large parameter space. 

In the previous section, we saw that the spatial scale 
of the expanding SN can affect significantly the emergent 
spectrum, mostly from the variation of the photospheric 
density. Let us assume that the model is fully ionized 
and that the continuum opacity is mostly due to electron- 
scattering (n e ~ 0.34 cm 2 /g), and investigate how model 
properties are affected by a change in the characteristics of 
the density distribution. Along a radial ray, the continuum 
optical-depth is given by r(r) = J r n e p(r')dr' so that at 

the photosphere, we have r p hot = K "^°l" 1 - ( R ^ j 

Let us take another model identical in all parameters 
except for those that characterize the density distribution 
(n 1 , po, n ', -Rphot, n') an d compare these at the photosphere: 



obtain 



Po.nR* ( R* 



1 \ -Rpho 



P0,; 



'7L 



R* 



1 \ -Rpho 



Within this simple approach we can determine how the 
photospheric radius changes if we keep the base density 
(and i?*) the same for both models, i.e. po,n = po,n'- We 



R 



phot.n' 

R* 



n-1 

n' - 1 



R 



phot,n 

R* 



Thus, for identical base density, if n' > n, Rphot.n' < 
fiphot,n, i-e. the photosphere moves in to hotter and slower 
regions. 

We show in Fig.^Jthe synthetic spectra of models dif- 
fering solely in density exponent n - in particular the base 
density is kept constant - covering values of 6, 8, 10 and 
12. We clearly see that as the exponent is increased, the 
spectrum shows a more and more pronounced UV flux, 
metal lines of decreasing strength in the optical but a 
strengthening Hei5875A. As expected, the photospheric 
radius decreases from 1.9, to 1.81, 1.54 and 1.4 times the 
base radius (kept constant) as one goes from n = 6 to 
12. Given the adopted Hubble law, a similar decrease in 
photospheric velocity follows (from 7,600 to 7,250, 6,150 
and 5,600kms _1 ). But there is a significant rise of the ef- 
fective temperature over that range, from 6,400, to 7,900, 
8,100 and 8,250K. These effects are reversed if instead one 
maintains the photospheric radius rather than the base 
density when changing the density exponent (Eastman & 
Kirshner 1989). 

Besides the general influence on the spectral energy 
distribution, changing the density exponent has also a sig- 
nificant effect on lines. For example, there is a systematic 
blue-shift of the observed Ha P-Cygni profile in Type II 
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Fig. 11. Illustration of the influence of the density exponent n on the synthetic flux distribution for n — 6, 8, 10 and 
12. All other model parameters are kept identical. The model atom used is small, including only the most abundant 
species, i.e. Hi, Hei, Cn, Ni, Oi, Fen and Fein. 
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Fig. 12. Illustration of the influence of the density exponent n on the synthetic flux distribution for n = 6,8, 10 in 
a cool model where hydrogen recombines near the base. All other model parameters are kept identical. Note that the 
most striking effect is on the set of Ca II lines near 8500A, the line appearing broader and stronger for smaller density 
exponents [color]. 
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SN, which can be modulated in models by changing the 
density exponent (Sect.[5J. Here we illustrate the influence 
the density exponent has on the strength and width of the 
Can feature near 8500A. This set of lines is present only 
in models that are cool enough, i.e. those where hydro- 
gen recombines at the photosphere. We show in Fig. IT21 
the synthetic flux distribution of three models differing 
only in density exponent (n = 6,8 and 10). Contrary to 
the previous figure, the effect on the overall spectrum is 
small, mostly because in this parameter space, the line 
and continuum formation is essentially contained within 
the (localized) region of the recombination front. However, 
as shown in Fig. ITS! the Can lines at 8498, 8542 and 8662A 
form over the entire outflow and are therefore extremely 
sensitive to the density gradient even above the recom- 
bination front, resulting in a broader and stronger Can 
feature for smaller exponents. However, Ha is more con- 
fined to the photosphere region and is thus only weakly 
sensitive to the outer density distribution. This explains 
why this set of Call lines can be used to reveal the pres- 
ence of outflow asymmetry, as inferred, e.g., by Kasen et 
al. (2003) in the Type la SN2001el, from the presence 
of strong polarization (with a different angle) in this Ca n 
spectral region compared to the rest of the spectrum. This 
line can therefore be used as a diagnostic for the density 
distribution and/or departure from sphericity of SN out- 
flows in general, when ionization conditions in the outflow 
are suitable, i.e. for effective temperatures at and below 
8000K. 

4.6. Lines of selected species 
4.6.1. Hydrogen lines 

At present there is an open debate concerning the capabil- 
ity of model atmospheres to fit the hydrogen Balmcr lines 
of Type II SN during the photospheric phase. Some groups 
fail to reproduce these lines (Schmutz et al. 1990 and 
Mitchell et al. 2002, both for SN1987A) while others fit 
them very well (Eastman & Kirshner 1989; Hoflich 1988, 
again for SN1987A). Schmutz et al. invoked the possibil- 
ity of clumping although they did not investigate whether 
it could in reality cure the problem, while a decade later, 
Mitchell et al. invoked energy deposition following 56 Ni 
mixing in the outer envelope of SN1987A, as soon as 4 days 
after core-collapse. Those who fit the hydrogen Balmer 
lines well for the same dates do not use clumping and 
ignore energy deposition. 

This situation is unsatisfactory since it is difficult to 
establish whether the problem is physical or numerical. 
We have shown in the previous section that CMFGEN re- 
produces extremely well the hydrogen Balmer lines, even 
though we overplot the synthetic spectrum and adopt a 
linear flux scale, two elements that have the potential 
to make fit discrepancies very clear. Although we show 
in Sect. [3] only three cases, all synthetic spectra com- 
puted for SN1987A and SN1999em fit the observations 
very well, provided that a strong hydrogen recombination 
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Fig. 13. Illustration of the relative line emission occur- 
ring at a given height in the outflow, both for one of the 
three components of the Can feature at ca. 8500A, i.e. 
Can8662A (solid line) and Ha (broken line). We use a 
color coding to differentiate between models with different 
density exponent. Note how extended the line formation 
region of the Can line is compared to that of Ha. The 
other two components of the Can triplet show a similar 
behavior and are thus not shown. £ is defined such that 
line flux is proportional to J £g? log r [color] . 

front does not develop, i.e. as long as hydrogen remains 
mostly ionized in the outflow. Indeed, if such a front ap- 
peared, CMFGEN always predicted Ha too weak by a 
factor of 2-3 compared to observations, both for the emis- 
sion and absorption components (see for example the Ha 
region in the cool models shown in Fig. I12fl . 

In CMFGEN, whose principal purpose is to determine 
the atmospheric structure and level populations, we sim- 
ply adopt a Gaussian shape, of constant width, to describe 
the intrinsic profile. In practice, the width of the profiles 
is set by the adopted micro-turbulent velocity. 

To save computational time, we have chosen a turbu- 
lent velocity wturb of 100 km s -1 . For fully ionized models, 
varying this value by a factor 2-3 has no effect on the 
emergent spectrum. However, for cooler models, we find 
that reducing this turbulence to 20kms~ 1 is important in 
order to reproduce better the profile shape and strength of 
Ha. In those cool models, a lower turbulent velocity gives a 
SED with a higher flux in the UV, and essentially a higher 
apparent outflow ionization (in fully ionized models, the 
effect on the SED is unnoticeable) . We understand this as 
resulting from the reduced line-blanketing, itself stemming 
from the reduced spectral band-width over which metal 
lines block the ionizing flux emanating from the thermal- 
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isation layer at the base. Similarly, late-stage simulations 
performed with fewer metal species show an appreciable 
increase in the Ha strength, again likely to result from the 
modulated amount of line-blanketing. 

4.6.2. Helium lines 

We have seen in Sect. 14.31 (Fig. 19")) that varying the metal- 
licity, and hence the magnitude of the line-blanketing, 
changes the strength of the He 1 5875A line. This He I 
line has proven to be difficult to model, even impos- 
sible without invoking unreasonable enrichments in he- 
lium (Schmutz et al. 1990, Eastman & Kirshner 1989). 
However, as can be seen from the detailed computations of 
Mitchell et al. (2002) for SN1987A and Baron et al. (2000) 
for SN1999em, PHOENIX reproduces very well this fea- 
ture in the early spectra taken. This seems to suggest that 
treating in non-LTE the metal lines that are responsible 
for the line-blocking is important. 

We show in Fig. [H CMFGEN fits to Hei5875A ob- 
served in SN1999em on the 30th of October and on the 
1st, 5th and 8th of November 1999 (Leonard et al. 2002a). 
The quality of the fits is very good, with only a modest 
helium enrichment (H/He = 5 by number). 

4.6.3. Nitrogen lines 

Based on a direct analysis with SYNOW, Baron et al. 
(2000) proposed that the feature blue- ward of Hei 5875A 
in the spectrum of SN1999em taken on the 30th of October 
was due to Nil. However, they could not confirm this 
statement with the more detailed analysis done with 
PHOENIX. They also highlighted the presence of a kink in 
the blue wing of H/3, which they do not associate with N II, 
but rather to some peculiar emission process occurring in 
the outflow. Similar intriguing features were observed in 
the early spectra of another Type II SN, SN1999gi by 
Leonard et al. (2002b). They observed the presence of un- 
identified features in the blue wing of both He 1 5875A and 
H/3, at blue-shifted velocities of ca. 20,000 to 30,000 km s" 1 
although, these have now clearly-defined P-Cygni profile 
shapes. Given the relatively higher density exponent char- 
acterizing the outflows of Type II SN at earlier dates, it 
seems unlikely that much emission could arise at veloc- 
ities (or distances), a few times the value at the photo- 
sphere, especially in recombination lines since these have 
a density-square dependent emissivity. 

We show in Fig. ^| synthetic fits to observations of 
SN1999em taken on the 30th of October 1999 (Leonard 
et al. 2002a). The dark-blue curve corresponds to the 
emergent spectrum when all species are included. It fits 
well the Hei5875A and H/3 regions, although a slightly 
higher expansion velocity would improve the fit to the 
He I line (i.e., a shift of the line formation region to larger 
velocities). Note the presence of kinks at around 4600A 
and 5500A, i.e. in the blueward region of each of these 
two lines. We over-plot the emergent spectra with bound- 
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Fig. 14. Synthetic fits (broken line) to observations (solid 
line) of SN1999em taken on the 30th of October (A), 
and the 1st (B), 5th (C) and 8th (D) of November 1999 
(Leonard et al. 2002a), showing the He 1 5875A line region. 
Note how well CMFGEN manages to reproduce the ob- 
served profile using only a modest helium enrichment over 
the cosmic value (H/He=5 by number), in agreement with 
stellar evolutionary predictions for Blue/Red supergiants, 
progenitors of Type II SN. This capability of CMFGEN 
is in stark contrast with failures to reproduce this He I 
line in the works of Schmutz et al. (1990) and Eastman & 
Kirshner (1989), suggesting the importance of treating all 
species in non-LTE. 



bound transitions exclusively of hydrogen (turquoise), he- 
lium (green) and nitrogen (red curve), normalizing all 
curves to the observed flux at 8000A. Unlike PHOENIX, 
CMFGEN is capable of reproducing the Nil line in the 
blue wing of He 1 5875A, but together with this line, it also 
predicts the equally strong N 11 feature around 4600A. In 
our model, these correspond to multiplets of Nn(3d-3p) 
centered around 4480A and 5490A, as well as N n(3p-3s) 
at 4630A and 5680A. 

The reason for the different results between PHOENIX 
and CMFGEN, regarding the N 11 lines and the kink in the 
blue wing of H/3, is unclear. Our fits to hydrogen Balmcr 
line profiles are much better, in particular the width and 
depth of the troughs, giving us confidence that we are 
not predicting spurious absorption too far to the red of 
line centers. Moreover, their Fig. 6 suggests that not just 
H/3 shows a kink beyond the point where the absorption 
trough meets the continuum, but that this feature is also 
present at a similar location for Hei5875A and Ha (note 
that their use of a very thick line makes this difficult to 
see) . Such a kink is however not observed in Ha (see their 
Fig. 3). 
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Fig. 15. Zoom into the optical spectroscopic observations 
(black) of SN1999em taken on the 30th of October 1999 
(Leonard et al. 2002a) showing synthetic fits for a model 
comprising all species as given in Table 1 (blue) or only 
bound-bound transitions of hydrogen (turquoise) , helium 
(green) and nitrogen (red). We strongly support the idea 
that the observed features around 4600A and 5500A are 
due to N II lines [color] . 



Our model fits show unambiguously that nitrogen lines 
are present in the optical spectrum of the early photo- 
spheric phase of SN1999em. In SN1999gi, spectroscopic 
observations taken after the 12th of December 1999 are 
very similar to those of SN1999em, so that we reach the 
same conclusion. Note however that for the first spectrum 
of SN1999gi (taken on the 10th of December 1999), the 
feature at ca. 4600A is much stronger than that at ca. 
5500A, a situation we have not been able to reproduce. 
By the 12th, the feature has both receded and weakened 
and is then well reproduced (and identified as an N n line) 
by our model for SN1999em on the 30th of October 1999 
(Fig.ll5(l. In our model, we use a factor ten enrichment of 
nitrogen compared to cosmic while we find that these N II 
features become unnoticeable when this factor becomes 
smaller than around 3. This strongly supports the idea 
that the progenitors of SN1999em and SN1999gi had core- 
processed nitrogen at their surface at the time of collapse. 

Interestingly, although the first spectrum of SN1987A 
shows the high ionization conditions for the presence of 
these Nil features, these are not observed. Together with 
the fact that a modest helium enhancement permits a nice 
fit of the spectrum, this suggests that the pre-collapse 
progenitor of SN1987A had a surface abundance closer 
to the current LMC abundances. Hence, identification of 
N II lines is a significant clue for enhanced nitrogen abun- 
dance in the outer layers of the progenitor star. Note that 
the over-abundance of nitrogen should be accompanied by 
a significant depletion of carbon and oxygen if it resulted 
mostly from the peeling-off of the outer envelope through a 
stellar wind mass loss. In the case of mixing, which can for 
example be made more efficient through stellar rotation, 
a more modest depletion in carbon and oxygen would re- 
sult. Unfortunately, we have not been able to identify lines 
from either of these species in the photospheric phase of 



Type II SN, although the addition of Ci in our model 
atom could alter this conclusion. We defer until a future 
study a proper discussion on this issue. 

5. Line formation in expanding outflows of hot 
stars and Type II SN 

Before concluding on these premiminary results with 
CMFGEN, we present a study on the formation of P- 
Cygni line profiles in Type II SN outflows, which reveals 
subtle differences from that observed, e.g., in the epony- 
mous star. Indeed, a striking feature in early-time spectra 
of Type II SN is the large blue-shift, by as much as ca. 
100-150A (5000-7000 km s" 1 ), of the Ha line profile com- 
pared to the rest wavelength of 6562. 79A. At later times, 
this blue-shift is less obvious, partially because the pho- 
tospheric velocity quickly decreases, but it is nonetheless 
present (see, e.g., Pastorello 2004). 

To investigate the origin of this feature, we first focus 
on the formation of scattering (resonance) and recombi- 
nation lines in Wolf-Rayet outflows. Their density distri- 
bution is much flatter than that of SN, following roughly 
a 1/r 2 law (imposed by mass conservation); non-LTE ef- 
fects are important but line-formation occurs over much 
larger spatial scales compared to the photospheric radius. 
Therefore, these objects offer a nice counterpart to line- 
formation conditions in SN, which we are trying to under- 
stand. 

We are particularly interested in determining where 
in the outflow the emergent flux originates from. Using 
models for nitrogen-rich Wolf-Rayet (WN) stars similar 
to those studied elsewhere (e.g., Hillier & Millier 1998, 
Herald et al. 2001), we show in Figs. llfiflTI grav-scale im- 
ages in the (x,p) plane of the flux- like quantity p ■ I(p), 
where x = (A/Ao — 1) c, A is the wavelength, p is the im- 
pact parameter and I{p) is the specific intensity (at A, or 
x). Note here that p is in units of the hydrostatic radius 
R* of the WN star under consideration. In Fig. ^| (WN5 
case), we show the case of the resonance line C IV 1548 A 
in a hot hydrogen-free WN model with an asymptotic ve- 
locity Woo = 2000kms _1 , while in Fig.lT7lfWN8 case), we 
show the case of the recombination line Ha in a cooler 
(hydrogen-rich) WN model (with Voo — 840 km s -1 ). In 
each figure, the top panel shows the synthetic line profile, 
which, for any scaled velocity x, corresponds to the inte- 
gral over all p of the quantity p ■ I(p) plotted in the lower 
panel. Thus, for any x, we directly identify the contribu- 
tion at p of the relative fraction made to the total line 
profile flux. 

The gray-scale illustrations contain a lot of informa- 
tion. First, beyond the velocity limits of the line (limited 
naturally to ztuoo), we primarily observe continuum pho- 
tons, although some line photons are observed outside this 
range (primarily on the red side) due to the effects of inco- 
herent electron scattering. In Fig. 1161 the continuum pho- 
tons originate at p values up to about pn m = 5i?», which 
is the location of the photosphere in this WN star model. 
In Fig. E| we observe a similar effect, although the pho- 
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tosphere is much closer in, at ca. 2i?*. But in both cases, 
the line remains optically-thick out to much larger values 
of p - a few tens of the corresponding . This is consis- 
tent with the generally-accepted fact that line formation 
in Wolf-Rayet outflows occurs over large scales compared 
to the photospheric radius. 



X n =1548 A 




Fig. 16. Hot star model (WN5). Bottom: Grayscale 
image of the quantity p ■ I{p) as a function of p and scaled 
wavelength x = (A/Ao — l)c, where p is the impact parame- 
ter (in units of the hydrostatic radius R*), I{p) the specific 
intensity along p (at x). Here, Ao corresponds to the rest 
wavelength of C IV 1548 A and c is the speed of light. The 
terminal velocity of the model is 2000 km s -1 . Top: Line 
profile flux, directly obtained by summing p- I(p) over the 
range of p. Thus, the line flux at x in the top panel corre- 
sponds to the cumulative sum of all contributions p ■ I(p) 
at x shown in the bottom panel, giving a vivid illustration 
of the sites at the origin of the observed line profile. 

Although the lower panels of Figs.ll6landll7llook qual- 
itatively similar, Fig. 1161 shows a P-Cygni profile, while 
Fig. 1 171 shows an emission profile. Close inspection of the 
lower panels shows that emission beyond p\i m is strong, 
roughly symmetric about line center (modulo turbulence 
effects): A key difference is that emission extends up to 
Wfrac ~ 0.5«oo for the WN5 model (Fig. El but right up 
to Vfc ac ~ Vqo for the WN8 model (Fig. ITTjl . In both cases 
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Fig. 17. Hot star case (WN8) Same set as in Fig. [TJl 
but for the Ha line in a WN8 model with a terminal ve- 
locity of 840 km s -1 . 

there is a strong flux deficit on the blue side for p < pn m , 
which is compensated for in the line profile, but only up to 
i Wfrac, by the line flux emitted at larger impact parame- 
ters. For the WN5 case, absorption beyond — 1000 km s -1 
is not compensated by emission and a trough appears in 
the resulting profile. For the WN8 case, absorption merely 
introduces skewness in an otherwise pure emission profile. 

Let us now consider the case of a Type II SN, taking 
a model similar to that presented in Sect. 13.21 We show 
the equivalent of Figs. 1161171 for that SN case in Fig. [TBI 
for Ha. i?* now refers to the base radius of the simula- 
tion also called i?o- The difference with the Wolf-Rayet 
case is large. Line and continua are confined to similar re- 
gions, both becoming optically-thin beyond pn m — 2i?». 
Interestingly, both emission and absorption occur below 
Pii m - Emission thus only appears on the blue side of the 
profile, which explains why the resulting emission profile 
is strongly blue-shifted. This occurs because the density 
drop-off is high in Type II SN (n = 10 in this model) 
so that emission is favored in the inner denser regions, 
strengthening the impact of disk-occultation. 

What controls the different behaviors shown in 
Figs. 1161171181 is the radial variation of the line source 
function Si compared to the continuum source function 
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Fig. 18. Type II SN model Same set as in Figs. 1161 
1171 for the Ha line formed in a SN outflow with similar 
properties to those described in Sect. 13.21 

at a radial optical depth of approximately r c = 2/3. We 
show in Fig. EH plots of Si/S c (t c = 2/3) versus contin- 
uum optical depth and versus radial velocity for the three 
models. Consider a photospheric (core) ray. The contin- 
uum intensity emanating from that ray is approximately 
S c (t c = 1) (w S c (t c = 2/3)), where r c now refers to the 
optical depth along the ray and we have made use of the 
Eddington-Barbier relation. On the other hand, the inten- 
sity in an optically thick line at some frequency in the line 
will simply be the line source function at the resonance 
zone for that frequency. For the WN5 star we see that the 
line source function is less than the reference continuum 
source function for r c < 2/3 and V > 900 km s _1 , and 
drops to very small values for V > 1200 km s _1 . Thus we 
expect no emission (as compared to the continuum) for 
V > 1200 km s _1 , as indeed seen in Fig.UHl For the WN8 
star the Ha source function is less than S c (t c — 2/3) be- 
yond 300 km s -1 , but it drops, by at most a factor of 4 
at 750 km s _1 . Thus the flux is reduced relative to the 
continuum, but never goes black. The reduced flux is eas- 
ily compensated by emission from rays not striking the 
photosphere, and hence we see no P-Cygni profile. 

The SN model is different yet again. From Fig. El we 
see the line source function drops below the continuum 
source function for velocities greater than 7000 km s . 
Thus, for such velocities, we expect to see reduced emis- 
sion, as seen from Fig. El Further, because the continuum 



opacity is mainly due to electron scattering, we can see 
emission arising from optical depths greater than t c = 1 
(V < 7000 km s _1 ). Because the line source function is 
larger than the continuum in this region, the line is in 
emission. This is again seen in Fig. ITSl Unlike the two 
WR models, we get significant emission for rays striking 
the photosphere. Because (by definition) the photosphere 
is thick we do not see much emission from the backside 
of the flow. Consequently the emission, for photospheric 
rays, is strongly skewed to the blue (Fig. 1 1 Hp . 




Fig. 19. Radial and velocity variation of the logarithm 
of the line source function normalized to the continuum 
source function at a continuum optical depth of two third 
(see text for details). We show the cases corresponding to 
Figs. 1161171151 i.e. for Civ 1548 in the WN5 model (top), 
Ha in the WN8 model (middle) and Ha in the SN model 
(bottom). 



Finally, we show in Figs. 1201 and 1211 the cases for SN 
models that have different density exponents, i.e. n = 8 
and n — 14 respectively compared to the model shown 
in Fig. El Indeed, the higher n is, the faster the density 
drops and the more blue-shifted the resulting Ha profile 
looks. While the source function follows a similar trend 
in these two additional cases (not shown), it is the modu- 
lated influence of disk occultation that gives rise to differ- 
ent blue-shift magnitudes in the resulting profile. Overall, 
the large density exponent characterizing the early evolu- 
tion of Type II SN is also a key ingredient that precludes 
the possibility of observing a clear departure from spheric- 
ity (if at all present), the line and continuum formation 
regions overlapping in space and being so confined. 
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Fig. 20. Same as Fig. [^J for a model with n = 8. 




Fig. 21. Same as Fig^|but for a model with n = 14. 



6. Conclusion 

We have presented first results based on the model at- 
mosphere code CMFGEN (Hillier & Miller 1998) for the 
spectroscopic modeling of photosphcric-phase Type II su- 
pernovae. To facilitate the modeling, several changes were 
made to the code: these include the ability to use a Hubble 
velocity law and a power-law density distribution, an 
adaptive grid to handle the steep H recombination front in 
cooler SN models, and an improved algorithm to compute 
the gray temperature distribution in a rapidly expanding 
outflow. Still to be included in our treatment are rela- 
tivistic effects, energy deposition from radioactive decay, 
and the possibility to base our computation on a com- 
puted hydrodynamical structure. We demonstrated the 
ability of CMFGEN to reproduce a large number of fea- 
tures in high S /N multi- wavelength spectroscopic observa- 
tions of Type II SN, using the very well-observed SN1987A 
and SN1999em. The observed spectral energy distribution 
from UV to the near-IR (when available) can be matched 
very accurately: the fit quality is in general within the 
10% level, although it is difficult to match observations 
in the UV range to better than 20-30%. We show exam- 
ples taken over a wide range of outflow ionization proper- 
ties, covering fully- (early-time) and partially- (late-time) 
ionized spectra. This good behavior of CMFGEN gives 
us hope that a wealth of information can be retrieved 
through a detailed analysis of multi-epoch observations 
of Type II SN. Based on a representative set of observa- 
tions, we show that Hei lines in early-time spectra can 
be well fitted by invoking modest helium enrichment typ- 
ical of an OB-supergiant progenitor. The non-LTE radia- 
tive transfer solution also predicts the presence of lines in 
the optical from higher ionization than before expected, 
in particular Fein and Nil. The latter provides a simple 
explanation for the appearance of "bumps" in the blue- 
wing of H/3 and He 1 5875A, generally interpreted in terms 
of over-dense emitting (in the corresponding line) struc- 
tures lying far above the photosphere. Overall, this study 
shows that taking into account as accurately as possible 
the enormous effects of line-blanketing is key for a de- 
tailed modeling of Type II SN. We also discuss through a 
number of different stand-points the problem of line and 
continuum formation in Type II SN, giving vivid illustra- 
tions for both. In particular, for the first time, we give 
an illustration to explain the sometimes strong blue-ward 
wavelength shift of, e.g., Ha in Type II SN. The shift 
arises from the combined effects of disk-occultation and a 
line formation region confined to the photosphere. 

Forthcoming studies will present more focused anal- 
yses on specific issues. Since most investigations done on 
Type II SN are with the aim of using them for the distance 
calibration of the Universe (as holds for la's), we will next 
study the Expanding Photosphere Method and the analy- 
sis of Eastman et al. (1996), following their approach but 
using CMFGEN to search for possible differences, in, e.g., 
correction factors. We will also discuss methods for deter- 
mining the photosphcric velocity, an issue of importance 
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since the rate of expansion of the radiating photosphere 
is one of the key quantities entering any such distance 
determination. Long term goals are to design a very reli- 
able method to determine distances in the Universe based 
on Type II SN, with application to a sample of very well 
observed objects. Analyses of individual objects is also 
sought, in order to constrain the evolutionary status of 
the SN progenitor based on a spectroscopic approach, to 
complement inferences made exclusively from photomet- 
ric measurements performed on the pre-explosion object 
(see, e.g., Smarrt et al. 2004). 
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